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Abstract.  A comparative  experimental  study  of  large-scale  space-time  structures  embedded  into 
the  broad-band  multi-scale  turbulent  plasma  fluctuations  of  the  purely  toroidally  magnetized  exper- 
iment TEDD1  is  presented.  Two  different  plasma  sources  arc  compared:  Plasma  can  be  produced 
cither  by  a thermionic  discharge  or  by  inductive  heating  by  a helicon  antenna.  Turbulent  plasma 
potential  and  plasma  density  fluctuations  were  analyzed  by  conditional  averaging.  It  turned  out  that 
large-scale  space-time  structures  develop  on  both  plasmas,  which  propagate  with  local  E x B drift. 

A contribution  of  diamagnetic  drifts  is  not  observed.  In  the  thermionic  discharge  dipole-like  density 
structures  are  found,  in  the  helicon  discharge  the  detected  structures  are  generally  of  monopole-like 
shape  and  are  considerably  sized  when  compared  to  the  thermionic  discharge.  Direct  measurements 
of  the  fluctuation  induced  transport  can  be  related  to  the  specific  properties  of  the  detected  large- 
scale  fluctuation  structures.  It  is  concluded  that  the  injection  of  the  negative  space  charge  via  the 
heated  filament  greatly  alters  the  properties  of  the  turbulent  fluctuations.  Although  the  basic  statisti- 
cal properties  are  substantially  unaffected  by  the  plasma  source  a strong  influence  on  the  formation 
of  large-scale  fluctuation  structures  and  the  associated  transport  is  found. 

INTRODUCTION 

It  is  a well  knows  feature  of  turbulent  systems  that  they  develop  large-scale  space-time 
fluctuation  structures.  Fluctuation  structures  can  be  observed  not  only  in  plasmas  but 
also  in  turbulence  in  neutral  fluids,  atmospheric  turbulence  etc.  In  plasmas  the  formation 
and  properties  of  space-time  fluctuation  structures  have  been  attracting  great  attention 
because  they  are  suggested  to  contribute  significantly  to  the  fluctuation  induced  transport 
perpendicular  to  the  magnetic  field.  The  large  fluctuation  amplitudes  can  lead  to  large 
transport  events  and  the  large  spatial  extent  of  the  structures  leads  to  transport  over  large 
radial  distances. 

Experimentally,  investigations  of  large-scale  fluctuation  structures  where  conducted 
in  linear  [1, 2,  3,  4]  as  well  as  in  toroidal  magnetic  field  geometry  [5,  6,  7]. 

In  linear  geometry  long-living  large-scale  monopole-like  structures  embedded  into  the 
turbulent  plasma  fluctuations  were  detected.  In  toroidal  geometry  large-scale  fluctuation 
structures  were  also  found  but  with  different  basic  properties.  Here,  the  detected  struc- 
tures have  dipole-like  shape.  The  interpretation  of  this  difference  is  complicated  because 
the  effect  of  the  magnetic  field  geometry  in  the  toroidal  devices  are  generally  strongly 
entangled  with  effects  arising  from  charge  sources  present  to  produce  the  plasma  in 
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FIGURE  1.  Schematic  drawing  of  the  TEDDI  device  (a)  together  with  the  two  plasma  sources,  the 
filament  source  (b)  and  the  helicon  source  (c). 


these  campaigns.  Based  on  numerical  simulations  it  was  suggested  that  the  formation  of 
dipole-like  structures  can  be  mainly  attributed  to  the  presence  of  the  charge  source  [8,  9], 
To  separate  the  geometry  from  plasma  source  effects  an  experimental  study  conducted 
in  a device  with  simple  toroidal  magnetic  field  geometry  is  presented,  in  which  two  dif- 
ferent plasma  sources  are  used.  First,  plasma  can  be  produced  by  the  injection  of  highly 
energetic  electrons,  which  is  the  standard  method  for  plasma  production  in  such  devices. 
Alternatively,  plasma  can  be  produced  ’quasineutral’  by  inductive  heating  without  the 
injection  of  net  charge. 

This  paper  is  organized  as  follows:  After  the  introduction  of  the  experimental  setup 
the  diagnostics  and  data  analysis  method  of  conditional  averaging  arc  briefly  described. 
Subsequently  the  experimental  results  of  conditional  averaging  for  both  plasmas  and 
direct  measurements  of  the  fluctuation  induced  transport  arc  presented.  Finally,  the 
results  are  summarized  and  conclusions  arc  drawn. 


EXPERIMENTAL  SETUP 

The  experiments  were  conducted  in  the  toroidal  plasma  device  TEDDI  [10,  1 1,  1 2 J.  A 
schematic  drawing  of  the  device  and  the  plasma  sources  is  shown  in  Fig.  I.  It  consists 
of  a toroidal  vacuum  chamber  with  three  wedge-shaped  segments,  where  the  plasma 
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TABLE  1.  Device  and  typical  plasma  parameters  valid  for  the 
filament  and  helicon  operation. 


Parameter 

Value 

large  device  radius 

0.3  m 

small  device  radius 

0.1  m 

neutral  Argon  gas  pressure 

3.2  ■ !0~2  Pa 

magnetic  field 

0.15T 

plasma  density 

1 • 1016  m~3 

electron  temperature 

4 eV 

ion  temperature 

0.03  eV 

filament  operation: 

discharge  voltage 

100  V 

discharge  current 

0.2  A 

helicon  operation: 

rf  frequency 

28MHz. 

rf  power 

200W 

source  is  placed  and  measurements  of  poloidal  cross  sections  are  performed.  The  device 
is  magnetized  by  12  planar  poloidal  magnetic  field  coils.  Consequently,  the  magnetic 
field  is  purely  toroidal.  The  typical  magnetic  field  strength  on  the  device  axis  is  0.1 5T. 
Plasma  can  be  produced  by  two  completely  different  discharge  mechanisms:  by  direct 
injection  of  high  energetic  electrons  (filament  discharge)  or  by  inductive  coupling  (he- 
licon discharge).  For  the  filament  discharge  (Fig.  1(b))  a tungsten  wire  is  inserted  into 
the  vacuum  vessel.  The  wire  is  directly  heated  by  a DC  current  and  the  emitted  electrons 
are  accelerated  between  the  negatively  biased  filament  and  the  grounded  poloidal  lim- 
iters placed  on  both  sides  of  the  source  region.  Plasma  is  produced  by  ionizing  collisions 
between  the  primary  electrons  and  the  Argon  neutral  gas  background. 

Alternatively,  the  plasma  can  be  produced  inductively.  For  this  discharge  mode  a 
helicon  antenna  [13]  (double-half-turn  right-helical  antenna)  is  installed  at  the  position 
of  the  filament  (Fig.  1(c)).  Measurements  of  the  magnetic  wave  field  suggest  that  a 
helicon  m — I mode  is  driven  by  the  antenna  at  a rf-frequency  of  10-30MHz  [10]. 
With  a rf-power  of  200W  similar  plasma  density  and  electron  temperature  are  obtained 
when  compared  to  the  filament  discharge.  1 The  basic  device  and  plasma  parameters  are 
compiled  in  Tab.  1.  The  important  difference  between  the  two  discharge  types  is  that 
a strong  negative  source  of  charge  is  present  in  the  filament  generated  plasma  whereas 
plasma  production  is  ’quasineutral’  in  the  helicon  discharge.  This  difference  facilitates 
for  direct  investigation  of  the  influence  of  net  sources  of  charge  on  the  formation  and 
properties  of  large-scale  fluctuation  structures. 


1 It  should  be  noted  that  helicon  discharge  in  this  context  does  not  indicate  the  very  efficient  antenna- 
plasma  coupling  that  leads  to  100%  ionization  degree  [14] 
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DIAGNOSTICS  AND  DATA  ANALYSIS 


Measurements  are  performed  with  electric  probes  exclusively.  Equilibrium  plasma  pa- 
rameter profiles  are  obtained  by  evaluation  of  Langmuir  probe  characteristics.  For  den- 
sity fluctuation  measurements  Langmuir  probes  are  negatively  biased  and  the  ion  satura- 
tion current  fluctuations  are  taken  as  proportional  to  density  fluctuations.  For  plasma  po- 
tential fluctuations  emissive  probes  are  used,  which  are  fairly  insensitive  against  hot  elec- 
tron components  [15,  16].  Fluctuation  induced  transport  measurements  are  performed 
with  triple  probes. 

Large-scale  fluctuation  structures  are  extracted  out  of  the  noise-like  turbulent  fluc- 
tuations by  the  statistical  method  of  conditional  averaging  (CA).  In  its  simplest  form 
CA  is  based  on  two  simultaneously  recorded  time  series:  On  fluctuation  time  series  <t»rcl 
is  recorded  at  a fixed  spatial  point  ro  and  acts  as  the  reference  signal  for  the  analysis. 
Another  fluctuation  time  series  <t>pos  is  recorded  at  an  arbitrary  second  spatial  position 
ro  + Ar,  here  called  the  positioned  dignal.  By  a predefined  condition  p on  OrC|  charac- 
teristic events  in  the  time  series  are  detected.  Together  with  a time  interval  At  around 
these  time  instants  r,  an  ensemble  of  sub-time  scries  out  of  <l>mov  of  length  N are  col- 
lected, which  are  taken  as  statistical  independent  realizations  of  a random  process.  Con- 
sequently, ensemble  averaging  of  the  sub-time  series  extracts  the  part  of  fluctuations, 
which  is  consistently  repeated  in  the  sub-time  scries  (coherent  part)  whereas  the  inco- 
herent fluctuations  are  suppressed.  The  procedure  reads: 

1 n 

^conCA)  + Ar,  x)  = — [dv,s(m  T Ac,  t,  + x)|<hrci  (r(1,/,')  = p]  . (1) 

/v  /•=  I 

The  idea  of  conditional  average  is  as  follows:  If  a large-scale  fluctuation  structure  passes 
the  position  of  the  reference  probe,  it  produces  a specific  signal  pattern  in  <l>r0f,  which  is 
detected  by  the  condition  p.  If  the  structure,  depending  on  its  lifetime  and  propagation, 
reaches  the  positioned  probe,  a similar  signal  pattern  is  produced  in  <hpos.  By  taking  the 
time  instants  /,  as  the  reference  time  for  the  analysis,  the  phase  relation  of  structures 
occurring  in  the  reference  and  positioned  signal  is  conserved  and  ensemble  averaging 
extracts  the  structure  form  the  superimposed  random  turbulent  background  fluctuations. 
To  obtain  a confidence  level  for  CA  a conditional  deviation  is  defined  as  the  standard 
deviation  at  each  time  lag  Ax  of  the  difference  between  the  conditionally  averaged 
structure  and  the  single  realization  of  the  sub-time  series.  This  quantity  is  normalized  to 
the  standard  deviation  of  the  reference  signal 

°caM #M) ' 01 

A value  of  Oca  well  below  unity  indicates  that  the  single  realizations  are  reproducible 
whereas  values  close  to  one  indicates  that  the  single  realizations  strongly  deviate  from 
the  average  result. 

For  the  current  analysis  an  amplitude  condition  of  la  together  with  a rising  slope 
condition  is  used.  The  slope  condition  ensures  that  events  are  always  found  on  the  same 
slope  of  the  fluctuation  signal  and  artificial  phase  shifts  due  to  the  finite  widths  of  the 
fluctuations  are  avoided. 


684 


filament  discharge 


helicon  discharge 


0.06 

0.04 

0.02 


-4  -2  0 2 4 -4  -2  0 2 4 

density  [a]  density  [o] 


skewness=  0.16 
kurtosis  = 0.03 


FIGURE  2.  Density  fluctuations  (top  row)  with  the  respective  probability  distribution  function  (middle 
row)  and  power  spectral  densities  (bottom  row)  for  the  filament  (left  column)  and  the  helicon  discharge 
(right  column) 


EXPERIMENTAL  RESULTS 

Basic  statistical  properties  of  turbulent  fluctuations 

In  Fig.  2 time  series  of  plasma  density  fluctuations  are  shown  for  both  discharges 
together  with  their  probability  distribution  function  (PDF)  and  their  power  spectral 
density.  The  density  fluctuations  in  both  discharges  show  an  irregular  behavior.  In  the 
filament  discharge  fluctuation  are  slightly  burst-like,  whereas  in  the  helicon  discharge 
the  density  fluctuation  show  a noisy  behavior.  This  is  also  expressed  by  the  PDF,  which 
is  peaked  and  asymmetric  for  the  filament  discharge.  This  is  quantified  by  the  non-zero 
higher  moments  skewness  and  kurtosis  (note  that  a Gaussian  distribution  has  s=0,  k=0 
in  this  representation).  In  the  helicon  discharge  the  PDF  is  of  almost  perfect  Gaussian 
shape.  Both  power  spectral  densities  show  a power  law  decrease  S(f)  f~a  over  one 
frequency  decade,  which  indicates  that  both  plasma  are  in  a fully  developed  turbulent 
state  [17].  The  decrease  is  stronger  in  the  filament  plasma  with  a spectral  index  of 
a = 6.5  whereas  in  the  helicon  plasma  the  spectral  index  is  a = 2.9. 


Conditional  averaging  results 

In  Fig.  3 typical  results  of  CA  of  the  reference  and  positioned  probe  are  presented 
for  both  discharge  types.  The  reference  signals  show  a pronounced  positive  fluctuation 
structure,  which  is  mainly  induced  by  the  positive  amplitude  condition  met  at  x = 0.  The 
temporal  width  of  the  pattern  is  larger  for  the  filament  plasma  when  compared  to  the 
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FIGURE  3.  Conditional  averaging  signal  at  a single  spatial  point.  Shown  are  the  conditionally  aver- 
aged reference  signals  (a),  (d),  the  conditionally  averaged  positioned  signals  (b),  (e),  and  the  respective 
conditional  deviations  as  defined  in  Rq.  2 (c),  (f)  for  the  filament  and  the  helicon  discharge 


helicon  plasma.  Stronger  differences  are  found  in  the  positioned  signal.  Here,  a dipole- 
like behavior  is  observed  in  the  filament  plasma  with  negative  and  positive  amplitudes 
of  similar  magnitude.  In  the  helicon  plasma  the  conditionally  averaged  positioned  signal 
shows  a positive  amplitude  fluctuation  with  only  little  adjacent  minima.  The  conditional 
deviation  is  zero  at  the  time  instant  where  the  condition  is  fulfilled  and  increases  with 
increasing  absolute  value  of  the  time  lag  t.  It  is  noticeable  that  the  temporal  confidence 
interval  is  larger  in  the  filament  plasma  than  in  the  helicon  plasma. 

To  obtain  space-time  information  of  large-scale  fluctuation  structures  measurements 
of  density  fluctuations  are  performed  in  a two-dimensional  section  perpendicular  to  the 
magnetic  field  lines.  The  results  are  shown  in  Fig.  4 for  the  filament  generated  plasma 
and  in  Fig.  5 for  the  helicon  plasma. 

Filament  discharge.  In  Fig.  4 the  results  of  CA  of  density  fluctuations  for  a complete 
poloidal  section  are  shown  in  normalized  values.  The  time  interval  extends  over  160;/.?. 
The  cross  at  t = 0 indicates  the  position  of  the  reference  probe.  Large-scale  minimum 
fluctuation  structures  with  adjacent  maximum  structures  are  found  in  the  complete  time 
window  under  consideration.  Thus,  the  observed  structure  always  appear  in  a dipole- 
like  configuration.  The  minimum  structure  propagates  radially  outwards  during  the  first 
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FIGURE  4.  Conditional  averaging  results  in  a two-dimensional  plane  for  density  fluctuations  in  the 
filament  discharge.  The  ambient  magnetic  field  points  into  the  plane.  The  solid  white  circle  indicates  the 
inner  limiter  radius.  Minimum  and  maximum  structures  are  emphasized  by  white  and  black  contour  lines, 
respectively,  which  represent  a level  of  30%  of  the  peak  amplitudes.  The  labeling  of  the  radial  and  vertical 
axes  arc  in  units  of  the  inner  limiter  radius  a. 


half  of  the  time  window  t = -80... 0 /.is  and  poloidally  during  the  second  half  from 
t = 0. . . 80 /as.  The  same  propagation  properties  are  found  for  the  maximum  structure. 
Comparison  of  the  propagation  properties  with  the  equilibrium  potential  profile  yields 
that  the  propagation  of  the  structures  is  fully  determined  by  background  ExB  drift, 
which  is  well  distinguishable  from  diamagnetic  drifts  due  to  the  strongly  asymmetric 
plasma  potential  and  plasma  pressure  profiles  [10,  12].  The  lifetimes  of  the  detected 
structures  are  larger  than  the  chosen  time  interval  and  the  spatial  extent  are  of  the  order 
of  the  plasma  radius. 


Helicon  discharge.  A different  picture  is  obtained  in  the  helicon  discharge.  The  re- 
sults of  conditional  averaging  are  shown  in  Fig.  5 in  the  same  representation  as  Fig.  4. 
Here  a negative  density  fluctuation  structure  develops  at  x = 0 — 80/i.v,  propagates  radi- 
ally towards  the  device  center  and  decays  at  x = -20 /js.  At  the  same  time  a maximum 
structure  rapidly  grows  at  x — —20. . . 0 /.is,  is  turned  at  the  center  of  the  device  and  de- 
cays at  x = 40 . . . 60 /.is.  It  is  visible  that  both  structures  do  not  exist  simultaneously.  Thus, 
the  structure  are  of  monopole-like  shape.  The  lifetime  is  significantly  reduced  and  the 
structures  are  sized  when  compared  to  the  structures  found  in  the  filament  plasma. 
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FIGURE  5.  Conditional  averaging  results  of  density  fluctuations  for  the  helicon  plasma  in  the  same 
representation  as  Fig.  4 

Fluctuation  induced  transport 

An  important  parameter  for  the  fluctuation  induced  transport  associated  with  space- 
time  fluctuation  structures  is  the  phase  relation  between  plasma  potential  and  plasma 
density  fluctuations.  To  obtain  the  phase  information  for  the  large-scale  fluctuation  struc- 
tures under  investigation  the  CA  procedure  was  repeated  in  the  two-dimensional  spa- 
tial domain  with  plasma  potential  fluctuations  as  reference  signal.  Similar  structures 
were  found  as  shown  in  Figs.  4 and  5,  but  phase  shifted  by  « n/ 2,  which  results 
in  a strong  fluctuation  induced  transport  caused  by  the  structures  [12].  Fig.  6 shows 
the  two-dimensional  amplitude  coherence  between  plasma  potential  and  plasma  den- 
sity fluctuations,  the  directly  measured  fluctuating  transport  at  the  plasma  outside,  and 
its  probability  distribution  function.  The  amplitude  coherence  show  a broad  structure, 
which  strongly  suggests,  that  the  formation  of  structures  is  based  on  a MHD  instability 
[18].  This  is  consistent  with  the  observed  large  phase  shift  between  plasma  density  and 
plasma  potential  fluctuations.  The  fluctuating  transport  and  their  respective  PDF’s  for 
the  filament  and  helicon  plasma  show  significant  differences.  In  the  filament  generated 
plasma  the  transport  is  burst-like  with  large  transport  events  directed  outwards.  The  PDF 
is  peaked  and  asymmetric  towards  large  transport  events.  In  the  helicon  plasma  a more 
noise-like  transport  is  observed  with  a symmetric  PDF.  This  can  be  directly  related  to 
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FIGURE  6.  Two-dimensional  amplitude  probability  distribution  function  between  plasma  density  and 
plasma  potential  fluctuations  (top  row),  time  series  of  the  measured  fluctuation  induced  transport  (middle 
row)  and  the  respective  PDF  (bottom  row). 


the  dynamical  properties  of  the  observed  large-scale  structures.  In  the  filament  plasma 
the  spatial  region  of  existence  of  structures  is  mainly  outside  and  the  phase  relation  be- 
tween density  and  potential  structures  lead  to  a transport  caused  by  these  structures  that 
is  directed  outwards.  This  strongly  suggests  that  the  large  transport  events  are  associated 
with  these  structures.  In  the  helicon  plasma  structures  exist  in  the  plasma  center  and  the 
phase  shift  between  plasma  density  and  plasma  potential  fluctuations  is  such  that  the 
resulting  transport  is  directed  towards  the  device  center  [12].  Consequently,  no  exposed 
large  transport  events  are  found  and  the  PDF  is  symmetric. 


CONCLUSIONS 

The  presented  experimental  study  was  addressed  to  the  question,  what  determines  the 
specific  properties  of  large  scale  turbulent  fluctuation  structures.  By  operating  the  de- 
vice with  two  completely  different  plasma  sources  geometry  effects  could  be  separated 
from  source-induced  effects.  Large-scale  structures  were  found  in  both  plasmas,  which 
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propagate  with  local  background  ExB  drift.  In  contrast  to  linear  magnetic  field  ge- 
ometry, where  generally  small  phase  shifts  between  large-scale  plasma  potential  and 
plasma  density  fluctuations  arc  observed,  the  magnetic  curvature  in  toroidal  magnetic 
field  geometry  leads  for  both  plasma  sources  to  a large  phase  shift  of  tc/2.  This  is  con- 
sistent with  numerical  simulations,  which  includes  both,  curvature  driven  instabilities 
and  drift  waves,  and  measurements  of  the  amplitude  coherence  between  plasma  density 
and  plasma  potential  fluctuations.  But  the  plasma  source  has  a great  influence  on  the 
specific  shape  and  lifetime  of  structures.  If  negative  space  charge  is  injected  into  the 
plasma  by  the  filament,  structures  generally  appear  dipole-like  and  have  long  lifetime. 
In  the  helicon  discharge,  where  plasma  is  produced  quasincutral,  the  structures  are  of 
monopole-like  shape  and  their  lifetime  is  considerably  reduced  when  compared  to  the 
filament  generated  plasma.  Specific  features  of  the  measured  fluctuation  induced  trans- 
port can  be  attributed  to  the  large-scale  structures.  In  the  filament  plasma  the  propagation 
and  phase  shift  of  structures  lead  to  transport  events  that  dominate  the  transport  time  se- 
ries. In  the  helicon  plasma  no  significant  contribution  to  the  transport  radially  outwards 
due  to  the  detected  structures  is  expected  and  consistently  the  measured  transport  is  more 
noise-like. 

The  results  indicate  that  it  is  not  sufficient  to  include  only  the  magnetic  geometry  in 
the  interpretation  of  particular  properties  of  turbulent  plasma  fluctuations.  As  it  was 
explicitly  demonstrated  for  large-scale  fluctuation  structures,  the  plasma  source  can 
strongly  alter  the  turbulent  fluctuations  and  the  associated  transport. 
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